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Abstract. The limit behavior of trajectories of dissipative quadratic stochastic oper- 
■ ators on a finite-dimensional simplex is fully studied. It is shown that any dissipative 

quadratic stochastic operator has either unique or infinitely many fixed points. If dissi- 
pative quadratic stochastic operator has a unique point, it is proven that the operator 
' is regular at this fixed point. If it has infinitely many fixed points, then it is shown that 

uj — limit set of the trajectory is contained in the set of fixed points. 
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^} . 1. Introduction 

(N 

Many situations in population genetics can be described by discrete dynamical systems 
[1],[6],[7],[9]. More precisely, dynamical systems in population genetics are generated by 
nonlinear maps from finite-dimensional space into itself. S.N.Bernstein in [1] introduced 
a class of dynamical systems which is now called a quadratic stochastic operator (q.s.o. in 
short). Q.s.o. has the following form: 



{V%)k — ^ Pij,k%i%j j 

where x = (x\,X2, ...,x m ) S S m ~ l and 



in 



S™" 1 = {x G R m : Xi > 0, ^ = 1} 

i=l 

is m — 1 dimensional standard simplex. Here the coefficients Pij,ki called heredity coeffi- 
cients, satisfy the following conditions 

m 

Pij,k = Pji,k > 0, ^Pij,k = i. 
k=l 

If the above conditions are satisfied, then one can easily check that q.s.o. maps the 
simplex into itself. 

Research on quadratic stochastic operators was later developed by a number of authors 
(see for example [2], [5], [13], [15]). The main problem in this area is the study of limit 
behavior of the trajectory of an initial point taken from the simplex. Note that the 
limit behavior of trajectories of quadratic stochastic operator was fully classified on one 
dimensional simplex by Yu.I.Lyubich in [8], [9]. However, the problem is still open in 
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higher dimensions (even in two-dimensional simplex). In some papers, classes of q.s.o. 
were considered with additional properties that enables the study easier (see [12]). For 
example, in [2], the class of Volterra q.s.o. was outlined which satisfies the additional 
condition Pij k = whenever k ^ {i,j}- It was shown that the trajectory of any nonfixed 
initial point, approaches the bound of the simplex. 

A q.s.o. is called doubly stochastic if Vx -< x for all x G S 1 " 1-1 , where " -< " is the 
notation of majorization, that is the classical ordering by comparison of partial sums of 
the coordinates after the nonincreasing rearrangement. The class of doubly stochastic 
operators was introduced in [3], [4], and the structural properties of the set of doubly 
stochastic operators were studied. The class of doubly stochastic quadratic operators is 
sufficiently large, therefore the study of limit behavior is complicated. However, by altering 
Vx and x in the definition of doubly stochastic operator, we find the class of dissipative 
q.s.o. The class of dissipative q.s.o. was introduced in [14], in which it was studied some 
properties of dissipative q.s.o.. Here we note that the definition of dissipativity is not 
given in terms of wandering set as in [11], but in terms of majorization. 

From the definition of dissipative and doubly stochastic q.s.o. one can see that the 
intersection of the set of dissipative and doubly stochastic q.s.o. 's is nothing but permu- 
tation operators of the simplex. Besides, in [14], it was shown that only identity operator 
can simultaneously be Volterra and dissipative. So, we assert that the class of dissipative 
operators is the new class of q.s.o., which makes the study of this class more interesting. 

We say that a q.s.o. V satisfies the ergodic theorem ([16]), if for any x G S"™ -1 the limit 

s + ViH W n ~ x x 

hm 

n— >oo n 

exists. In [14] it was shown that dissipative operator satisfies the ergodic theorem. How- 
ever, the convergence of Cesaro mean of the trajectory of dissipative q.s.o. does not imply 
the convergence of the trajectory of the dissipative q.s.o. Therefore, the study of limit 
behavior of the trajectory of dissipative q.s.o. is still open. 

In the present paper, we classify the limit behavior of dissipative q.s.o. on finite dimen- 
sional simplex. The paper is organized as follows. In section 2 we give some terminology 
and notations. Also, we provide some previous results concerning dissipative q.s.o. In 
section 3 we give main results of the paper. 

2. Preliminaries 

In this section we give some definitions and state some previous results. Let 

m 

S™" 1 = {x G R m : Xi > 0, ^xi = l} 

i=i 

be a (m - 1)— dimensional simplex. Then the vectors 

e fc = (0,0,--- - - ,0), 

k 

(k = 1, m) are its vertices. For a C / = 1, 2, • • ■ , m, the set F a = {x G S™ 1 ' 1 : Xi = 0, i 
a} is called a face of the simplex. 

For x G S" 1 " 1 Let us put x± = (arjij, xp], ■ ■ ■ £[m])> where (xjij, arp], • • • X[ m ])- nonincreas- 
ing rearrangement of (x±, X2, • • • x m ), that is x^ > xp] > ■ ■ ■ > X[ m ] . 

Definition 1[10]. We say that x is majorized by y on S" 1 ^ 1 , and write x -< y(or y y x) 

if : 

k k 

J2 x \i] <J2v\i], k = l,m - 1. 
i=i i=i 
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It is easy to see that for any x G S m we have 

(-,-,■•• ,-) <x< (1,0,0,- ■■ ,0). 
mm m 

The following proposition can be found in [10]. 
Proposition 2.1. The following statements are mutually equivalent 

k k 

1) x <y, that is £ x [{ \ < £ V{i\, k = l,m - 1. 

i=l i=l 

2) x = Py for some doubly stochastic matrix P. 

3) Vector x belongs to the convex hull of the ml vectors obtained by permuting the 
components of the vector y. 

Definition 2. An operator V : S" 1 ^ 1 — > S m ~ l is called dissipative if 

Vx yx, Vx G S m ~ l . (1) 

The reason that the operator with above property is called dissipative is simple. Since 
Vx >- x, then according to the above proposition the point Vx does not belong to the 
convex hull of the points, obtained by permuting the components of x. The same can 
be said about points V 2 x and Vx, V 3 x and V 2 x, and so on. So, if we consider points 
x, Vx, V 2 x, ■ ■ ■ , then these points dissipate inside the simplex. 

It was stated in [14] dissipative linear operators are just permutations. Well, it is not 
generally so. Here the 

Observation. Let us consider the case when V is a linear dissipative operator, that 
is Vx = Ax, where A = (ciij)i j=TTn is an m x m matrix. Since Vx >- x, then by putting 
x = e, we have Ae^ >- ej. At the same time, it is easy to see that Aei -< ej. That is why 
(Aei)i = ( e i)li which means that only one component of the vector Aei is 1 and the others 
are 0. Hence, the operator can either be permutation operator or maps the simplex into 
its face which makes the study of the operator easy in either way. Here the example of 
linear dissipative operator which is not permutation: 

V : {x 1 ,x 2 , x 3 ) -> (x 1 + x 2 , x 3 , 0). 

However, if we consider the interior of the simplex, then linear dissipative operators are 
only permutations. Therefore, it is more interesting to study non-linear dissipative ones. 
In what follows, we consider quadratic operators. Recall [9] a stochastic operator V : 
gm-i _^ gm-i ^ g ca n ec [ a quadratic stochastic operator if it has a form 

m 

( Vx )k = Pij,kXiXj, (2) 

where x = (x\,x 2 , ...,x m ) G S m ~ 1 . Here the coefficients pij^ satisfy the following condi- 
tions 

m 

Pij,k = Pji,k > 0, ^2?ij,k = 1- (3) 
k=i 

More information concerning q.s.o. can be found in [9]. We note that linear operator on 
S" 1 ^ 1 can be considered as a quadratic operator since we can always multiply the operator 

m 

by E x ' = 1- But, since the study of the limit behavior trajectories of linear dissipative 
i=i 

operators are simple, here and henceforth we exclude this kind of operators from the set 
of dissipative q.s.o. 



4 DISSIPATIVE QUADRATIC STOCHASTIC OPERATORS 

Let us consider the following operator 

(Vx)l = x\ + x\ + X% + X\X2 + XiX 3 + X 2 X 3 , 

(Vx) 2 = xix 2 + xix 3 , 
(Vx) 3 = x 2 x 3 . 

Simple calculations show that (Vx)\ > max{x±, x 2 , x 3 } and (Vx)i + (Vx) 2 > max{x\ + 
x 2 ,x\ + X3}, (Vx)\ + (Vx) 3 > x 2 + X3. Therefore, above operator is dissipative quadratic 
stochastic operator. 

From the definitions 1 and 2, it is easy to see that if an operator is dissipative, then 
by changing its components it preserves its dissipativity. This gives more examples of 
dissipative q.s.o. 

Now, for the purpose of completeness we state two Lemmas from [14] and provide their 
short proofs since we are going to use them frequently. 

Given q.s.o. V we denote djj = (j>ij,i,Pij, 2 , ■ " Pij,m) Vi,j = l,m, where pij^ are the 
coefficients of q.s.o. V. One can see that aij G S" 71-1 , for all i,j = l,m. 

Lemma 2.2. Let V be a dissipative q.s.o. Then the following conditions hold 

{au)i = ei Vi = l,m. 

Proof. Due to dissipativity of F one has Fx >- x, Vx € S" 71-1 . Now by putting x = ei 
we get ej -< Ve*. On the other hand, we have >- x, Vx G S 171-1 . That is why (Vej)j, = 
( e i)^ = ei. Then the equality Fe^ = implies the assertion. □ 

The above lemma implies that any dissipative q.s.o. can be written as 

(Vx) k = ^2 x i + 2 ^2Pij,k x i x j k = l~m, (4) 

iSctfc i<j 

where 

m 

a k C I = {1, 2 ■ ■ ■ , m}, Oi n aj = 0, i ^ j, J a fc = /. (4') 

fe=i 

Lemma 2.3. Zei 6e o dissipative q.s.o. 

(i) // j G a fco , i/ien p i: ,- jfeo = (a^j > §, Mi = l,m. 

(ii) Ifm>3, then (a^-Vi = 0, Vfc > 3, Vz = l,m. 

Proof, (i). Let j G afc and x = (1 — A)ej + Ae^. Here, as before, ej,ej are the vertices of 
the simplex and A is sufficiently small positive number. It is easy to see that x^j = 1 — A 
and (Var)[ij = (Vx)k - Since Vx >- x, then x^j < (Vx)^, so 1 — A < (Vx)k or 

1_ A < (l-A) 2 + 2^ iifc0 A(l-A). 

The last inequality implies that Pij,k > \ ■ 

(ii). Denote Pijk* = maxpyj. One can see that (Vx)k* = (a>ij)\2\- Now from 

X[l] + x [2] < (Vx) {1] + (Vx) [2 ] 

we obtain 

1 < (1 - A) 2 + 2(py, fco +^, fc *)A(l - A). 

From this inequality we get Pij,k a +Pij,k* > 2X \i-x) = 2(1- \) = + - L This y ields 
Pij,fc + Pij,k* = 1 and (aij)[ k] = VA; > 3, Mi = T~m. □ 
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Now, let us recall some terminology. Let a; € S m ~ l . Then the set {x°, Vx°, V 2 x°, • • • } 
is called the trajectory starting at the point x°. The point x° satisfying Vx° = x° is 
called a fixed. The set of all fixed points of the q.s.o. V is denoted by Fix(V) A q.s.o. 
V : S m ~ l — > S*™" 1 is called a regular if the trajectory of any x £ S m ~ l converges to 
a unique fixed point. Let V be a q.s.o. Then the set co(x°) = f] \J {V n x } is called 

fc>0n>fe 

co-limit set of trajectory of initial point x° € S m ~ l . From the compactness of the simplex 
one can deduce that 7^ for all x° G S 1 ™ 1-1 . 

We provide statements that were proved in [14]. 

Let us consider the case a\ = I and a k = for k 7^ 1. Then the operator has the 
following form 



{Vx)i = E ^ 2 + 2 E Pij,iXiXj, 

i=l i<j 

(Vx) k = 2 E Pij,kXiXj, 2 < k <m. 

i<j 



(5) 



Theorem 2.4. ^4 g.s.o. given by (5) is regular. Its unique fixed point is e±. 

Now we turn to another case, namely let a\ = /\{/} «2 = {1} and a k = 0, Vfc > 3. 
Then operator (5) has the following form 

(Vx) 1 = E x i + 2 E 

1=1, j^/ J<j 



(yx) 2 = xf + 2 e Pij,2XiXj, 

i<j 

(Vx) k = 2 E Pij,kXiXj, 3<k<m. 

i<j 



(6) 



Theorem 2.5. If I ^ 2 then the operator (6) is regular and has a unique fixed point e\. 
If I = 2 then the operator (6) has infinitely many fixed points. Moreover, oj-limit set of 
trajectory of any initial point a; belongs to co{e\, e2}, here coA denotes the convex hull of 
a set A. 

As we have seen, in the above two theorems the limit behavior of the trajectory of the 
operator was studied in some particular cases. That is why, it is still important to study 
the limit behavior of the trajectory of the operator in general case that will be done in 
the next section. 

3. Main results 

In this section the main results are given. We show that any dissipative q.s.o. has 
either a unique or infinitely many fixed points. We also study the limit behavior of the 
trajectory of the dissipative q.s.o. 

Theorem 3.1. Any dissipative q.s.o. has either unique or infinitely many fixed points. 

Proof. Since the operator is continuous, and maps convex and compact set into itself, then 
according to Bohl-Brouwer theorem there exists a fixed point of the operator. 
Further, let us write down a dissipative operator in the following form 

( Vx )k =^2 x i + 2^2Pij,kXiXj k = T~m, 

i£a k i<j 

m 

where a k C / = {1, 2 ■ ■ ■ , m}, a« PI aj = 0, i 7^ j, (J a k = I. 

k=l 

We prove the statement of the theorem by considering several cases. 
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Case 1. There is no k such that k £ a k . That is, k £ a k for all k = l,m 
Case 2. There exists a unique ko such that ko £ ak - 

Case 3. There are numbers fci, fo, ■ ■ ■ j &n(2 < n < m — 1) such that fcj G = 1, n. 

Now we prove the statement for the Case 1. Let us assume that k ^ for all k = l,m. 
Then there exist k±,k2,--- ,k n (n < m), such that k{ € a^o..), where it is a permutation 
of the set fci,&2,--- j^n and 7r(z) 7^ i, Vz due to /c ^ WLOG we may assume that 
fcj = i, Vi = 1, n. 

Since any permutation is a composition of cycles then n = tt\ o 712 o ■ ■ ■ o n p , where in 
are cycles with the length not less than 2. 

Let us consider the case when p = 1. WLOG we may assume that 7r is a backward shift, 
that is, = i — 1 for i = 2, n and 7r(l) = n. Then the operator has the following form 

(Vx) k = xl +1 + Y xf + XI 2pij,kXiXj, l<k<n-l' 
«ea fc \{fe+i} i<j 

(Fx) n =x\+ Y x i + E 2 Pij,n x i x j, \ r>j\ 

i&a n \{l} i<j 

{Vx) k = Y x i + 2 Y Pij,kXiXj, n + 1 < k < m. 

i€a k i<j 

Let l<k<n(iik = n then we set n + 1 = 1). Then 

(Vx) k = X 2 k+1 + ^ X i +^2 2 Pij,k x i x j = 
i€a k \{k+l} i<j 

= Xk+1 + ^2 X "i +S ^Z 2 Pij,kXiXj + X 2 k+1 - X k+ i = 
i&a k \{k+l} i<j 

= Xk+1 + ^2 ^ + ^2 2 Pij,kXiXj - ^ x iX k+1 . 
i£a k \{k+l} i<j i=T~m,i^k+l 

Thus, the operator (7) can be rewritten as 

(Vx) k = x k+ i+ Y xf + Y 2 Pij,kXiXj - JY x i x k+i, l<k<n-l^ 

i£a k \{k+l} i<j i=l,m,i^k+l 

_ m 
(Vx) n = Xi + Y X 2 + Y 2 Pij,nXiXj - Y x iX!, 
iea n \{l} i<j i=2 

(Vx) k = Y x i + 2 Y Pij,kXiXj, n + 1 < k < m. 

i&a k i<j 

(8) 

Let us put 

L k = 2 Pij,kXiXj - ^2 x i x k+l- (9) 
i<3 i=T^n,i^k+l 

We show that Lk > 0, \/k = 1, n. Indeed, since fc + 1 € a&, VA; = 1, n (with n + 1 = 1) 
then according to Lemma 2.3, we have 2pik+i,k > 1 Vi = l,n, i 7^ + 1. Therefore, 

£fc = 2 Plk+l,kX±Xk+l + 2p2k+l,k x 2Xk+l + • • • + 2pkk+l,kXkXk+l + 2pk+2k+l,k x k+2Xk + ■ ■ ■ 

+2p mfc+ i i £ : 2 m Xfc + i+ ^ 2pi jifc XiXj-(xiXfc + iH hXfeX fc+ i+x fe+ 2Xfc+iH hx m x fc+ i) = 

= (2pifc+i,fc - l)xix fe+ i + (2p2fc+i,fc - l)x2Xfc+i H h (2p fefc+ i ife - l)x fc Xfc+i+ 
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+(2pfc+2fc+i,fc - l) x k+2 x k H h (2p mfe+ i ife - l)x m x fc+ i + ^2 2pij,kXiXj > 0. 

Further, since ir is a permutation with maximal length, then 

£ x? = (10) 

iGQfe 

for n+1 < k < m. Indeed, the set ct k for n+1 < k < m does not contain numbers 1, 2, • • • , n 
because these numbers belong only to one of the sets a±,a2, • • • ,a n . Moreover, k £ a k 
according to the assumption. If jo £ a k for jo > n + 1 j / k + 1 then there exist numbers 
ji,j2,--- ,3r{jr > n + 1, r > 1), such that j £ a k , ji £ a jo ,---j r £ aij x , k £ a jr . 
Let ir' be a permutation such that ir'(k) = jo, vr'(jo) = ji, • • ■ j 7r'(jV— l) = jVj tt'OV) = 
then a = 7r o 7r' is the permutation on {1, 2, • • • , m} with greater length than 7r, which is a 
contradiction. That is why, a k = for n + 1 < k < m. Hence, numbers k + 1, k + 2, • • • , m 
belong to one of the sets a\, ct2, ■ ■ ■ ,a n . 

Since L k > for all k = l,n, then by letting Vx = x and by summing up (8) from 
k = 1 to n we get 



E E ** 2 = o. 

fc=l 

m m n 

Due to (4'), it is clear that Y Y x i = Y x i ■ Moreover, from (10) we get Y Y x 



fc=liG«fc k=l 



Y x f- Therefore, 



k=i 



fe=l i€a fc \{fc+l} 



EE^-E^ 2 = E^-E^ 2 = E 



fe=l 



fc=l 



fc=l 



i=n+l 



which implies x n+ \ = x n+ 2 = ■ ■ ■ = x m = 0. 

In addition, we obtain x\ = X2 = ■ ■ ■ = x n which means that (^, ^, • • • , ^, 0, • • • ,0) is 
a unique fixed point. 

Now, let p > 2. There is no loss in generality in assuming p = 2, that is ir = ir\ o tt 2 . 
WLOG we may put 7Ti (i) = i — 1 for i = 2, ni with 7Ti(l) = n\ and 112(1) = i — 1 for 
i = ni + 2, n-2 with ^(ni + 1) = «2- Then the operator has the following form. 



{Vx) k = x 2 k+1 + Y x i + Y 2 Pij,k x i x j, 1 < fc < ni - 1 
«ea fc \{fe+i} i<j 

(T^x) ni = X]_ + ^ Xj + ^ 2pij^ ni XiXj , 
•iea ni \{i} »<j 

(Fx) fe = x^ +1 + J] X? + Y,2Pij,k x iXj, n 1 + l<k<n 2 -l 

i€a k \{k+l} i<j 



(Vx] 



ri2 — ^ni 



+1 + S x f + S 2 Pij,n2 X i x ji 



iGa n2 \{"-l+l} 



(Vx) k = Y x i + 2 Y Pij,k x i x j, n 2 + 1 < k < m 

i€a k i<j 



(11) 



which can be rewritten as 



8 
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(Vx) k = x k+ i + Y x 1 + T, 2 Pij,kXiXj - J}2 XiX k+1 , 1 < k < ni - 1 

i£a k \{k+l} i<j i=l,m,i^k+l 



m 



(Vx) m = xi + Y x i + Y 2pij, ni XiXj - Y XiXi, 

iea ni \{l} i<j i=2 

(Vx) k = x k +i+ Y x i + Y 2 Pij,kXiXj - Y XiX k+1 , n 1 + l<k<n 2 -l 
i£a k \{k+l} i<j i=l,m,i^k+l 

{Vx) n2 = X ni -i-i + Yl X^ + Yl 2 Pij,ri2XiXj — Yl x iX ni -\-i, 

i6o;„ 2 \{ni+l} i<j i=T~m,i^ni+l 

(Vx)k = Y x i + 2 Y Pij,kXiXj, n 2 + 1 < k < m. 

(12) 

If we put Vx = x and sum up (12) from k = 1 to k = n 2 , we obtain 

712 ™2 ™2 

J2(Vx) k = x k + S 1 + S 2 + ^ L fe> 
fc=i fc=i fc=i 

ni rt2 

where Si = £ £ x? (m + 1 = 1) and S 2 = Y E a;? (n 2 + 1 = n x + 1). 

fe=liea fc \{fc+l} fe=ni+l iGa fc \{fc+l} 
"2 "2 

Since > VA; = l,n 2 and £ (T^x)fc = Y x k , then we have 

k=l k=l 

ri2 ri2 ri2 r*2 "2 

X] Xfc = ^(^fc = ^2 x k + Si + S 2 + ^2L k >^2x k + St + S 2 

k=l k=l k=l k=l k=l 

which implies Si + S 2 = 0. As we showed the relation (10), we can also show that Y x 



2 

i 

i£a k 

m 

2 



for n 2 + 1 < k < m. Therefore, = Si + S 2 = Y x h an d hence Xj = Mi > n 2 + 1. 

fc=ri2+l 



r 2 

i£a k \{k+l} 



Now, since Y x i an d L k are positive, we let Vx = x, and from (12) we get x k > 



x k+ i for k = 1, n\ (x ni +i = xi), which implies x± = x 2 = ■ ■ ■ = x m = a. Consequently, we 
get x fc > Xfc+i for k = m + l,n 2 (x n2+ i = x ni +i), which implies x„ 1+ i = x ni+2 = ■■■ = 
x n2 = (3. So, the set of fixed points is 

{(«,«,••• ,a,/3,/3,--- ,/3,0,0,--- ,0) GS™" 1 ! n lQ + n 2 /3 = 1} 

Thus, there are infinitely many fixed points. 

Case 2. Assume that there exist a unique k\ G a^, and WLOG we assume that &i = 1 
There are two possibilities 

(a) . There is a permutation tt of the numbers k 2 , k^, ■ ■ ■ , k n (2 < k% < m, 1 < n < m— 1) 
such that fcj G awfcj- 

(b) . There is no any permutation it of the set such that hi G a^.). 

In case (a), WLOG we may assume that ki = i and ir(i) = i — 1 for i = 3, n and 7r(2) = n 
then the operator has the following form 
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(Vx)i =x\+ Y xj + Y 2pij,iXiXj, 

ieai\{l} i<3 

(Vx) k = x 2 k+1 + Y x i + Y 2 Pij,kXiXj, 2<k<n-\ 

i£a k \{k+2} i<j 

(Vx) n = x\ + Y x i + Y 2 Pij,nXiXj, 
iSa„\{2} i<j 



(13) 



(Vx) k 

which can be rewritten as 



Y x l + 2 Y Pij,kXiXj, n + 1 < k < m, 

i&a k i<j 



(Vx)i=xi+ Y xj + Y 2pij,iXiXj - Y XiX\, 

i£ai\{l} i<3 i=2 

(Vx) k = x k+ i + Y xj + Y 2pij,kXiXj - _Y XiX k+2 , 2 < k < n - 1 

iea k \{k+2} i<j i=l,m,i^k+2 

{Vx) n = X 2 + Y xj + Y 2pij,nXiXj - Y x iX 2 , 
«Ga„\{2} i<j i=l,m,i^2 

Y X 1 + 2Y Pij,kXiXj, n + 1 < k < m. 



{Vx) 



«6ct fe 



Kj 



(14) 

We use the similar arguments as in Case 1 and deduce that the set of fixed points of 
(14) is 

{(1 - (n- l)a, a, a, ■ ■ ■ ,a,0,0---0)| < a < -} 

n— 1 times 

Let us now assume that (b) holds. The operator has the following form 



(V^x)i =x\+ Y xf+Y 2pij,iXiXj, 

ieai\{l} i<3 

(Vx) k = Y xj + 2 Y Pij,kXiXj, 2 < k < m, 

i£a k i<j 

which can be rewritten as 

m 

(Vx)l=Xl+ Y xj + Y 2 Pij,lXiXj - Y XiXl, 
iea-l\{l} i<j i=2 

(Vx) k = Y xj + 2 Y Pij,kXiXj, 2 < k < m. 

i£a k i<j 

If (b) holds, then it is easy to show (as was done in (10)) that 

xj = 0, VA; = 2~^i. 



(15) 



(16) 



i£a k 



m 

Hence, Y x j = Y x j ■ So, if we set Vx = x in (16), then it follows that x 2 = X3 = 

i6Ql\{l} i=2 

• • • = x m = 0, which means that e\ is a unique fixed point for (16). 

Case 3. Let us suppose that there are numbers k±, k 2 , ■ ■ ■ , k n (2 < n < m — 1) such 
that ki G «fc.,Vi = l,n. It is sufficient to show for n = 2 and k\ = 1, k 2 = 2. 

Then the operator has the following form 
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(Vx) 1 =x\+ Y X i + H 2p ijA XiXj, 

ieai\{l} i<j 

(Vx) 2 = xl+ Y X l+T, 2Pij,2XiXj, 
iea 2 \{2} i<j 

(Vx) k = Y xj + 2 Y Pij,kXiXj, 3<k<m, 

i&a k i<j 



which can be rewritten as 



(17) 



(Vx)i=Xi+ Y x i + I] Ipij^XiXj - Yl XiXi, 
i€ai\{l} i<j i=2 

(Vx) 2 = X 2 + Y X i + J2 2Pij,2XiXj - _Y XiXi, 
i£a 2 \{2} i<j i=T~^, i^2 

{Vx) k = Y x i + 2 E Pij,kXiXj, 3 < k < m. 

i€a k i<j 



(18) 



Now, if there exist numbers Is, h, ■ ■ ■ , l n (3 < h < m, 1 < n < m — 2) such that 
k G a-irHi), then by putting ^ = i and 7r(z) = i — 1 with 7r(3) = n we find that the set of 
all fixed points is 

{(q,/3, 7, 7 ,^- ,7 ,0,0---0)| a + /? + (n - 2) 7 = 1}. 

n—2 times 

If there is no such numbers then by using arguments as (10) was shown, we can also 
show that 

iSa fe 



0, VA; = 3, m, 



which means that Y x \ + Y x \ = Y x 1- So, if we put Vx = x in (18), then it 

i6ai\{l} iea 2 \{2} i=3 

follows that X3 = X4 = • • • = x m = 0, which means that the vectors Ae! + (1 — A)e2 are 
fixed points of (18). □ 

Remark 3.2. It can be seen from the proof of the theorem that the set of fixed points of 
operator has shown exactly up to permutation of coordinates. In those cases when operator 
has a unique fixed point, this fixed point is either the vertex of a simplex or the center of 
a face of the simplex. In those cases when operator has infinitely many fixed points, then 
the set of fixed points is some face of the simplex. 

Now let us consider some examples. 

Example 1. Consider the following operator on 2D simplex. 



(Vx)i = X2X3, 

(Vx) 2 =x\ + x\x 2 + X1X3, 

(Vx) 3 = x\ + X% + X\X 2 + X1X3 + x 2 x 3 . 

By putting Vx = x and solving the system of equations we find that (0, 0, 1) is a unique 
fixed point. 

Example 2. Consider the following operator on 3D simplex. 
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(Vx)i = x\ + X\%2 + X1X3 + (IX1X4 + X2X4, 
(Vx) 2 = X% + X1X3 + X 2 X 3 + 6X3X4, 

{ Vx) 3 = x 2 2 + x 1 x 2 + x 2 x 3 , 

(Vx)/i = (2 — a)xiX4 + (2 — 6)x3X4, 

where 1 < a, 6 < 2. The set of fixed points of the above operator is (1 — 2A,A, A,0), 
where < A < \. 

Example 3. Consider the following operator on 2D simplex. 

(Vx)i = x\ + x\ + X X X 2 + X1X3 + x 2 x 3 , 

(Vx) 2 = x 2 x 3 , 
(Vx) 3 = x\ + xix 2 + X1X3. 
Above operator has a unique point (|,0, \). 

In general, for any given vertex of the simplex or any given center of a face of the 
simplex, it is possible to provide an example of dissipative q.s.o. with a unique fixed point 
at the given point. 

Theorem 3.3. Let V be dissipative q.s.o. Then one of the following statements always 
holds. 

• The operator is regular. Its unique point is either a vertex of the simplex or the 
center of its face. 

• The operator has infinitely many fixed points, u— limit set of any initial point is 
contained in the set of fixed points i.e. uj{x) C FixiV). 

Proof. The way of the arguments in the proof is similar to the proof of the previous 
theorem. Here we divide the proof into 3 cases as well. 

Case 1. In this case the operator has one of the forms (8) and (11). Suppose that the 
operator is of form (8). Define 

tp(x) = xi + x 2 H h x n . 

Then 

XL XL Yl TL 

^x) = j>x), = j> + £ Yl ^ 2 + E L - ( 19 ) 

i=l i=l k=l iea k \{k+l} i=l 

Here L« satisfy (9). Since Lj > 0, then from (19) it follows that tp(Vx) > <p(x). Hence, 
the sequence {ip(V k {x)) : k = 1, 2, • • • } is non-decreasing and bounded. That is why the 
limit lim (p(V k x) exists. Let us put C = lim (p(V k x). 

k~ ¥00 k— >oo 

Taking (10) into account one can rewrite (19) as 

n n m n 

<p(vx) = y,( Vx \ = !>+ E ^ 2 + E L - ( 2 °) 

i=l i=l i=n+l i=l 

From (20) we get 

n n m n 

^{v k+i x) = Y,(v k+1 x) t = E(^ fe *)* + E ( yhx ^ + E L ^' = 

i=l i=l i=n+l i=l 

m n 

= ^{v k x)+ y, (v fc x)?+E L *' ( 21 ) 

i=n+l i=l 
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where 

L k = 2p ij , s (V k x) i (V k x) j - (V k x)i(V k x) s+1 s = l~n, k = 1, 2, ■ ■ ■ . 

1 <J i=l,m,iy^s+l 

m 

Since lim ip(V k+1 x) = lim ip(V k x) = C and L k > 0, then (21) implies (V kx )i ~> 
as k — > oo, Therefore, (V k x)i — s- for all i = n + 1, m, as k —> oo. Taking into account 

m _ 

^(F fe x)i = 1 we get C = 1. Furthermore, if k —> oo then ^ (^ fea; )? for 

«=i «ea s \{s+l} 

s = l,n and — >■ 0. Therefore, from 

(V k+1 x) l = (V k x) t+1 + (V k x)f + L^ l<i<n, 

igQi\{j+l} 

we get 

lim (V k x) i+1 = lim (V k x)i, 

k— ¥oo k— >oo 

which implies lim (V k x)i = ^ Vi = l,n. This means that the operator is regular with 

k^oo 

the unique fixed point (£, ^, • • • , 0, 0, • • • , 0). 

Now, let us assume that the operator is of form (11). Using the above method, we 
define 

lf(x) = X\ + X2 H h x n2 

and show that lim ip(V k x) = 1. Moreover, it can also be shown that 

k~ >oo 

lim (V k x) i+ i = lim {V k x)i 

k— 5>oo k^oo 

for i = l,m and 

lim (U fc a;)i + i = lim (F fc x)j 

fc— >-oo A;— >oo 

for i = ni + l,n2- The last equalities imply that the trajectory of any initial point ap- 
proaches the set of fixed points of the operator. 

Now, we turn to Case 2. Assume that (a) holds. Then we define 

(p(x) = x\ + x 2 H h x n 

and show that lim ip(V k x) = 1. Its proof was done in Case 1, so we do not repeat it. 

fc— s-oo 

Moreover, we have 

lim (V k x) i+ i = lim (V k x)i 

k— s-oo fc— s-oo 

for i = 2, n — 1. That is why, for any initial point x° G S" 71-1 we have 

w(x°) C {(a, /?,-•• 0, 0)| a + (n - l)/3 = 1} = Fix(V). 
Now, if (b) holds, then we put ip(x) = x\ and show that lim (p(V k x) = 1, which means 

k^oo 

that V is regular at e\. 

Let us consider the last Case 3 in which V is of form (18). If there do not exist numbers 
Z3, Z4, • • • ,l n (3 < li < m, 1 < n < m — 2) such that G a -K(k) f° r some permutation 
7T, then we define ^(x) = xi + x 2 and show that lim (p(V k x) = 1, which means that 

fc— s-oo 

w(x) C Fix(V). 
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Otherwise we put k = i and = i — 1 with 7r(3) = n. Then we define 

ip(x) = x\ + x 2 H hi„ 

and show that lim 92(V fc x) = 1. Furthermore, 

fc— s-oo 

lim (V k x) i+ i = lim (V fc x)j 

holds for i = 3, n — 1. Hence, 

w(x )c{(a,/3,7," - ,7,0,- •• ,0)| a + + (n - 2) 7 = 1} = Fix(F). 

□ 

Corollary 3.4. Let V be a dissipative q.s.o. on 2D simplex. 

If V is regular, then its unique fixed point is one of the following points: 

ei,e 2 ,e 3 , ^(ei + e 2 ), ^(ei + e 3 ), i(e 2 + e 3 ). 

// V has infinitely many fixed points, then either ui(x) C {Aei + (1 — A)e 2 | < A < 1} 
or uj{x) C {Aei + (1 - A)e 3 | < A < 1} or u{x) C {Ae 2 + (1 - A)e 3 | < A < 1}. 

4. Conclusion 

The main achievement of this paper is the classification of the limit behavior of the 
trajectory of dissipative q.s.o. on finite dimensional simplex. The main results are The- 
orem 3.1, Theorem 3.3 and Corollary 3.4. Most of the papers on quadratic stochastic 
operators are considered in small dimensions of the simplex. But, in the present paper it 
was considered in any finite dimensional simplex which made the investigation more dif- 
ficult. The methods, which were used for proving results can be used in other disciplines 
of mathematics, namely nonlinear analysis, dynamical systems and ergodic theory. Here 
some directions in which research can be carried on. 

1. One can consider polynomial or any continuous function V, defined on the simplex 
satisfying Vx >- x and study the limit behavior of its trajectory. Research can be started 
when V is a cubic operator, that is V(Xx) = X 3 Vx, A > 0. 

2. Dissipative quadratic stochastic operators can be defined on infinite dimensional 

oo 

space. Infinite dimensional simplex can be defined on i\ as S = {x G l\ : Xi > 0, Yl x i = 

i=i 

1}. Majorization for infinite cases is developed. But, the problem is still difficult since 
infinite dimensional simplex is not compact. We also can not state that dissipative q.s.o. 
has at least one fixed point. 
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